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The Choh Uhlenbeck equation was the first successful generalization of the Boltz- 
mann equation to include triple collisions in a systematic manner. Detailed anal- 
ysis of its triple collision term has led to the important concept of ring collisions, 
which create dynamic correlations between particles with a common collision his- 
tory. This review shows how the ring collisions have been at the root of most major 
developments in kinetic theory of the last 40 years, both for continuous fluids as 
well as for lattice gas automata. In present times they have become the standard 
tool to analyze long spatial and temporal correlations in simple and complex fluids. 



1 Historical perspective 

In the sixties many efforts were made to generalize the Boltzmann equation to 
dense liquids. The historical importance of the Bogoliubov Choh Uhlenbeck 
theory and the Choh Uhlenbeck equation is well expressed by the following 
quote from Brush' bookid on the history of kinetic theory: "At Michigan, Soon- 
Tahk Choh, a Korean student of Uhlenbeck, was doing apparently similar but 
ultimately more fruitful work. Early in 1958 he completed his dissertation, 
The kinetic theory of phenomena in dense gasesn, a work never published but 
frequently cited by later physicists who heard about it" . 

In an effort to put Choh's work in perspective I briefly review some im- 
portant events, relevant for the development of the kinetic theory of dilute and 
dense fluids. In a small variation on Cohen's classification in his Fifty years of 
kinetic theorya I divide the history of kinetic theory into four periods: 

• Classic Era 1855-1945 

• Renaissance 1946-1964 

• Modern Era 1965-1985 

• Post-modern Era > 1985 

The selection of important developments is based on Brush's boc 
overview in Choh's thesis, and on review papers of personal preference t 

The Classical Era started with Maxwell's derivation in 1859 of the ve- 
locity distribution in thermal equilibrium: the birth of kinetic theory. Soon 
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afterwards, in 1872 Boltzmann formulated his famous equation for the time 
evolution of the single particle distribution function, f{x,t) = /(r, v, i), in a 
dilute gas outside equilibrium. 



The collision term on the right hand side, which is quadatic in /, represents 
the change in / through binary collisions. It is based on the assumption of 
molecular chaos, expressing the absence of precollision correlations between 
two particles involved in a binary collision. In other words, the Boltzmann 
equation takes only sequences of uncorrelated binary collisions into account. 
This implies that the pair distribution function /'■-'^"^{xiX2,t), just before par- 
ticles 1 and 2 enter each others interaction range, factorizes into 



The Boltzmann equation shows the arrow of time; it is an irreversible equation, 
which is lacking the time reversal invariance of the mechanical equations of 
motion. In fact, f{x,t)dx represents the average number of particles in a 
small element, dx — drdv of phase space. The Boltzmann equation neglects 
fluctuations and can be considered as a mean field theory. Boltzmann was able 
to prove the i?-theorem and discuss the approach of f{x, t) towards thermal 
equilibrium, described by Maxwell's velocity distribution function. 

The remaining problems were to derive the Navier Stokes equation and to 
calculate the transport coefficients. The possibility of an explicit calculation of 
viscosity and heat conductivity was realized in 1917 when Chapman and En- 
skog developed the Chapman Enskog method for solving the Boltzmann equa- 
tion. This multi-time scale method assumes that transport processes, which 
occur on time scales large compared to the mean free time ioj are described 
by the normal solution, /(r,v,t) — /(v|n, u, T), which depends on space and 
time variables (r, t) only through the first five moments of the distribution 
function, a(r,t) = {n(r, i), u(r, i), T(r, i)}, i.e. the local density n, the local 
flow velocity u and the local temperature T . 

In this time regime, the temporal and spatial variations over a mean free 
time to and a mean free path — vt^, provide a small parameter, the non- 
uniformity^arameter, fi ^oV, which can be used as a systematic expansion 
parameters. Time derivatives, dtf ^ {df/da)dta, can be calculated by carry- 
ing out a //-expansion of the local conservation laws for a(r, t), where dta is at 
least of 0(/i). Then, the Boltzmann equation can be solved by a /z-expansion, 
i.e. / = /o -)- /t/i -I- . . where /o = n(m/27rA:T)''/2 exp[-TO|v - u\'^/2kT] is the 
local equilibrium distribution with d the dimensionality of the system. The 



atZ + v-v/ = /(//). 



(1) 



/(2)(a;iX2,t) ^ f{xi,t)f{x2,t) ~ /(r,v,i)/(r,V2,i). 



(2) 
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Chapman Enskog method yields a hnear integral equation for the deviation 
/i, which is proportional to the gradients of T and u, and provides explicit 
expressions for the heat conductivity and viscosity. 

In 1922 Enskog proposed a rather accurate '0, though approximate, gen- 
eralization of the Boltzmann equation, to account for transport processes in 
a dense fluid of hard spheres. This was achieved by replacing Boltzmann's 
molecular chaos assumption (2) by 

/(2)(xia:2,i) ^,g(|ri2| =a)/(ri,vi,i)/(ri-^,V2,i). (3) 

The local equilibrium function for two spheres at contact, ^(cr), accounts for 
static correlations in a hard sphere fluid, and the different space dependence of 
the two /'s accounts for the difference in position of the two colliding spheres, 
and represents the new transport mechanism of collisional transfer, i.e. in- 
stantaneous transfer of momentum and energy over the hard sphere diameter 
a through the interaction. This is in general the dominant transport mech- 
anism at high densities. The original Boltzmann equation takes only kinetic 
transport into account. 



2 BBGKY-hierarchy 



The Renaissance began in 1946 vdth Bogoliubov's article Problems of a dynam- 
ical theory, in statistical physicsu. The fundamental basis was the BBGKY- 



hierarchy la, which is named after Bogoliubov, Born, Green, Kirkwood and 



Yvon, and reads 



(at + vi-vo/i = Jdx2e,2f[l\ 

{dt+£{12))fil^ - /dX3(013+023)/l23 (4) 



and so on, where 6*12 is the weak scattering operator, 

dV{ri2) I ( d d 



712 = 



(5) 



dvi2 m \9vi (?V2 
with pair potential V{r), and two-particle Liouville operator 
£(12) = £"(12) - 012 = VI ■ Vi + V2 ■ V2 - Oi2- (6) 

"'rather accurate', that is in moderately dense systems where the mean free path £0 is 
larger than the hard core diameter a. In dense liquids, where £0 << a, the Enskog equation is 
inadequate for describing spatial variations with wavelength A in the range £0 < X < cr. These 
wavelengths are relevant for analyzing dynamic structure factors as measured in neutron 
scattering experiments. However, the revised Enskog equation of section 12 is very adequate. 
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This set of equations is an open hierarchy of coupled equations that expresses 
the time rate of change of the s-particle distribution function /'^''^ in terms of 
the (s + l)-particle distribution function, provided the interactions are pairwise 
additive. 

BogoHubov and Kickwood were interested in non-equihbrium properties of 
fluids and dense gases EJ. Uhlenbeck considered Bogoliubov as the originator 
of modern methods in kinetic theory, as he introduced the basic concepts of 
systematically approximating the dynamics of moderately dense gases in terms 
of small groups of two, three, four,... isolated particles, just like in Mayer's 
cluster expansion for the equation of state in thermal equilibrium. This set 
the stage in 1954 for Choh's thesis work. 

The main goal of this paper is to analyze the impact of the Bogoliubov 
Choh Uhlenbeck theory on later developments. Consequently I will refer to 
links with later developments in kinetic theory for Lattice Gas Automata ^ 
(LGA), in which the BBGKY-hierarchy has been used as wellllJ. 

An extremely important parallel development occurred also in 1954, which 
led to the Green Kubo formulas or time correlation function method for trans- 
port coefficients in fluids EJ. The shear viscosity, for instance, is given by 

^ = I dt{J,yit)Jxvm, (7) 



VksT Jo 

where Jxy{t) is the microscopic A^-particle momentum flux, which contains 
both kinetic as well as collisional transfer contributions. The generic form of 
a Green Kubo formula is a time integral over an equilibrium time correlation 
function, such as the stress correlation hmction in (7). It involves in fact 
two limits: first the thermodynamic limit nf / . .) /V is taken, followed by a 
time integral over an infinite time intervallllE3. Similar Green Kubo formulas 
for LGA were derived by Dufty and Ernst 113, in which the time integrals are 
replaced by time sums and the initial term in the summand has a weight one 
half. 



LGA are fully discretized models for non-equilibrium fluids Q, defined on a regular 
space lattice, where time is discrete, where point particles live on the nodes of a lattice, 
and have a small set of allowed velocities corresponding to the b nearest neighbor lattice 
vectors. The dynamics consists of a collision step, followed by a propagation step. If £ 
particles {£ = 2,3, ... ,6) are on the same node (zero range interactions), then they suffer 
£-tuple collisions according to well defined collision rules, which are obeying momentum 
conservation in fluid models, and lacking it in purely diffusive or reactive models. As the 
interactions in LGA are not pairwise additive, but involve genuine ^-tuple collisions, the right 
hand side of the first hierarchy equation for LGA involves already all distribution functions 

with £ = 2,3, . . . ,b + 1. Ofcourse, a similar structure occurs in continuous fluids with 
non-additive 6-body interactions. 
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3 Bogoliubov Choh Uhlenbeck theory 

This section will only give an outline of the basic ideas of this theory. The 
interested reader is referred to Cohen's lecturesEI about this subject. 

The system of interest is a moderately dense gas under conditions of 
standard temperature and pressure, composed of particles interacting through 
short range, strictly repulsive forces. Here the typical duration of a collision 
is To = 10""'^^ seconds, whereas the mean free time is to — 10~^ seconds. The 
new concept is that the existence of well separated time scales brings about a 
contraction of the description in terms of fewer variables, as time progresses. 

There exist three time stages: initial stage, kinetic stage and hydrody- 
namic stage. In the initial stage {t < tq) the full Liouville equation is required 
to describe the fast variation of the A^-particle distribution function. In the 
subsequent kinetic stage (tq < t < to), which is the first coarse grained time 
scale, the single particle distribution function f{x, t) is the only slow variable, 
and the higher order function Z*-"^ are assumed to be time independent Junc- 
tionals Z*^''-' (. . . \f{t)) of f{x, t), which implies a loss of memory about the initial 
state |. Once f'^'^\x,X2\f(t)) is determined, one has a closed kinetic equation, 
the generalized Boltzmann equation. The pcrturbative calculation proceeds 
by performing a formal density expansion. In the final hydrodynamic stage 
(t > to), which is the second coarse grained time scale, a further contraction 
of the description occurs. Here the only slow variables are the hydrodynamic 
fields a(r,t) = {n(r, t), u(r, t), r(r, t)} and f{x,t) becomes a time- and space- 
independent functional of these fields, f{r,t) = /(v|a(r, i)). This assumption 
is a generalization of the normal solution in the Chapman Enskog method. 
Once /(v|a(r, t)) has been determined, one has obtained a closed set of hydro- 
dynamic equations, which contains the Navier Stokes equations, and general- 
izations thereof, as well as explicit expressions for the transport coefficients. 
The pcrturbative calculation proceeds by performing a formal /i-expansion, 
where fi ~ is the non-uniformity parameter. 

How can these ideas be implemented to actually derive the generalized 
Boltzmann equation? Clearly, in order to solve an open hierarchy, like (3), 
some kind of closure relation is required. In the present theory closure is 
obtained by imposing the factorization assumption as a boundary condition 
on the BBGKY-hierarchy at t = — oo, expressing the absence of correlations 
in the infinite past. This is technically expressed through the relation 

lim S^r{l2)f^^HxiX2\S°_J{t)) = 

r — >oo 

"^As will be discussed later on, it is the assumption about the existence of a time- 
independent functional, to general order in the density, that breaks down and leads to the 
divergence of the density expansion of transport coefRcients, to be discussed in section 5. 
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x,(-t) x,(-t) 

Figure 1: Action of streaming operators. 

5(12)/(,xi,<)/(x2,t) = f{xl,t)f{x;,t), (8) 

and similar relations for higher order distribution functions. Here St {12) is 
the streaming operator for two interacting particles, St{l2) = exp[tC{12)], 
and 3^(12) is the free streaming operator, 5'5'(12) = cxp[t£o(12)], where the 
Liouville operators C and Cq are defined in (6). The streaming operators S't( 12) 
and 5°(12) generate the trajectories of two isolated particles with or without 
interaction respectively. Moreover <S(12) is defined through 

5(12) = lim 5_^(12)5°(12). (9) 

r— >oo 

Both streaming operators <S(12) and S'_t(12) act directly on the phases Xi 
with i = 1,2, by replacing them respectively by x* = S{12)xi and Xi{—t) = 
S-t{12)xi. The latter operator generates backward dynamics. The action of 
both streaming operators on xi and X2 is sketched in Figure 1. 

It is exactly the factorized boundary condition at t = —oo, which intro- 
duces the arrow of time into the Boltzmann equation and its generalizations. 
The BBGKY-hierarchy equations themselves - without the boundary condi- 
tion ~ are invariant under time reversal, as required by classical mechanics. 
The final ingredient in obtaining the generalized Boltzmann equation is the 
use of the small parameter density, n = N/V, as a systematic expansion pa- 
rameter. It enables one to perform an n-expansion of the unknown functionals 
/''*'' = /o*'' + "^./i"'' + • • ., of the equations of the BBGKY-hierarchy, and of the 
boundary condition (8). 
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4 Choh Uhlenbeck equation 



The implementation of these ideas were realized in Choh's thesis, and led to 
the Choh Uhlenbeck equation in 1958: 

{dt+v-V)f = I{ff) + Kifff) + ... (10) 
where I{f f) is the binary collision term, 

/(//)= J dx2e^2S{12)f{xut)f{x2,t), (11) 

and K{f f f) is the triple collision term, which will be specified below. It is an 
important new result, that generalizes the Boltzmann equation to include triple 
collisions. Bogoliubov had already shown in his 1945 article that /(//) for the 
case of hard spheres reduces to the collision term in the Enskog Boltzmann 
equation with .g(ri2| = a) = 1. This binary collision term I{ff) does take 
into account the difference in position of the colliding particles, the so-called 
coUisional transfer. In the limit of low density, it gives a higher order density 
correction to the Boltzmann equation. 

In the notation of Choh's thesis, the triple collision operator has the form, 

K{fff) = j dx2ei2 drS.rm J dxs[...]S°{12)Ul^J{x,,t) (12) 

with 

[•••] = y"d.T3[(0i3 + 023)5(123) -5(12)0135(13) -5(12)0235(23)]. (13) 

A more transparent form was obtained in the later cluster expansion methods 
of Green and Cohen (see review □), namely 

Kifff) = J dx2ei2 J dx3{Sil23) - Sil2)Sil3) 

- 5(12)5(23) +5(12)}/(xi)/(a;2)/(a;3). (14) 

It is instructive to analyze the different types of contributions from 5(123), 
illustrated in Figure 2. They contain genuine triple collisions, where three 
particles are simultaneously in each others interaction sphere. Once they have 
collided, the particles will separate and never recollide again due to the short 
range repulsive potential. In hard sphere systems, where interactions are in- 
stantaneous, this term is vanishingly small. 
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Figure 2: Genuine triple collisions (a) and ring collisions: recoUision (b), cyclic (c) and 

hypothetical (d) collision. 



Moreover, 5(123) contains sequences of two uncorrelated successive binary 
collisions (12) (13) and (12) (23). Such collisions have already been accounted 
for in the Boltzmann collision term /(//). So, they are subtracted out from 
the possible sequences of binary collisions among three particles in (14). 

Finally 5(123) contains sequences of three correlated successive binary col- 
lisions, such as the recoUision (12)(13)(12) and the cyclic collision (12)(23)(13), 
generally referred to as ring collisions. A more systematic analysis of the sep- 
arate terms may be found in Refu. 

In the second part of his thesis, Choh obtained the solution to the Choh 
Uhlenbeck equation in the hydrodynamic stage by following Bogoliubov's gen- 
eralization of the Chapman Enskog method, and derived formal expressions for 
the first density correction to the viscosity and heat conductivity for a mod- 
erately dense gas in three dimensions with short range repulsive interactions. 
The viscosity, for instance, has the form 

V = Vq + Vin + V2n'^ + . . . , (15) 

where 770 is the low density Boltzmann result and rji is the Choh Uhlenbeck 
contribution. It is expressed in terms of matrix elements of the triple collision 
operator, and contains contributions from genuine triple collisions, as well as 
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from the correlated ring collisions, illustrated in Figure 2. For the special case 
of hard spheres, Sengers and coworkers have computed 771, and found ±hat rji 
differs only a few percent from the predictions of the Enskog equationd. 

Similar ^calculations for the viscosity in LGA || have been carried out by 
Brito et al.llil. Here it was found that the ring collisions give corrections on 
the order of 10 % to the Boltzmann mean field results for all densities, and the 
resulting viscosities are in excellent agreement with those of computer simula- 
tions. I want to stress that the Boltzmann mean field equation and the ring 
kinetic equation for LGA are not restricted to low densities. The Boltzmann 
equation for LGA does include apart from binary collisions genuine triple, 
quadruple, • • • , 6-tuple collisions, as well as static correlations, accounting for 
the exclusion principle, which mimics the hard core interactions. 

The seven years after the appearance of Choh's thesis in 1958 may be de.- 
scribed as a period of consolidation. Using the cluster expansion method □, 
Green and Cohen derived the formal structure of the £-tuple collision terms 
{i = 2,3,4,.. .) in the generalized Boltzmann equation, in close analogy to the 
Mayer cluster expansion for reduced distribution functions in thermal equilib- 
rium, and it was believed that transport coefhcients could be calculated in the 
form of a virial expansion where the coefficients rje in (15) were determined by 
the dynamics of {£ + 2) isolated particles. Similar cluster expansion methods 
were applied to the Green Kubo formulas for viscosity and heat conductiv- 
ity. After some initial confusion about a possible difference between the Green 
Kubo formulas and the Choh Uhlenbeck equation, McLennan was the first to 
show that both methods yield the saine results for the first density correction 
to the transport coefficients in fiuidsEj. 

5 Discovery of the divergence 

The cluster expansion method of Green and Cohen or the alternative inversion 
method of ZwanzigEj made it feasible to investigate the fundamental assump- 
tion of the Bogoliubov Choh Uhlenbeck theory about the rapid decay of initial 
correlations within a time to- This assumption provided the formal justifica- 
tion for the rapid approach of the pair- and higher order distribution functions 
F^'^\xiX2 ■ ■ ■ X3, t) to the time independent functionals Z*^*-* {xiX2 ■ ■ ■ X:j\f{t)). 

In 1965 Dorfman and Cohen, Weinstock, as well as Goldman and Frieman 
were the first to make phase space estimates of the triple and higher order 
collision terms u, and to show that the virial expansion contains coefficients 

No reviews on the kinetic theory beyond the Boltzmann mean field equation are available 
for LGA. Therefore I refer in each section to the original LGA publications to show the 
parallel developments. 
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that are divergent with mcrcasing time; more specificany 



{ 



Vo + »7i(i)n- + ■ • ■ 



(d=2) 



(16) 



where the quadruple coUision contribution ?/2(i) in three dimensions and the 
triple collision ^contribution rii(t) in two dimensions are logarithmically diver- 
gent as t — > ooB. Kawasaki and Oppenheim were the first to sum the ring dia- 
grams. Subsequently the divergences and resummations were consolidated by 
more detailed calculations of van Leeuwen and Weijland, Swenson atLd-McLen- 
nan, Haines together with Dorfman and Ernst, Sengers, and FujitaHfl. 

It is of interest to elaborate on this fundamental discovery, which marks the 
advent of the Modern Era in kinetic theory, where collective effects determine 
the behavior of transport coefficients, even at low densities. To do so, I present 
the derivation of the log t divergence, as given in Ref. □ for hard sphere models 
and Lorentz gases. 

Consider the cca-integral in (14) at fixed xi,X2- Let J* dxs denote the 
phase space where the integrand of the recoUision term (12)(13)(12) is non- 
vanishing, and let a be the range of the repulsive interaction. The configuration 
with dx2, = dv^dr^ for the occurrence of a recoUision can be simply estimated 
in the coordinate system used in Figure 3. 



The solid angle dv^, in which particle 3 has to hit particle 2, such that 
particle 2 will recoUide with particle 1 after a time ti, is proportional to dva ~ 
{a / Vi2ti)'^~^ . For particle 3 to perform this well-aimed collision its position 
has to be inside a small collision cylinder dr^ ~ a'^'^v^dti. The total phase 




V. 




Figure 3: Phase space for recollisions. 
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space for the recollision event is then obtained by integrating over all possible 
intermediate times, yielding for large times t: 

This estimate illustrates that the triple collision term in the Choh Uhlenbeck 
equation in three-dimensional systems approaches a finite limit, implying that 
rji in (16) exists, and that the initial state is indeed forgotten, in agreement with 
Bogoliubov's assumption about the existence of a time- independent functional. 
However all three conclusions, valid for three dimensions, break down in two 
dimensions, where the Choh Uhlenbeck term diverges as logt. 

Next, I consider a similar phase space estimate for the ring collisions in 
the quadruple collision term, for instance (12) (23) (24) (12). The additional X4- 
integration yields an extra factor ti in (17). Consequently the relevant phase 
space, 

f d.. ^ £ dt.tA-^ ^ { ^ ^3-. 1^ ; I (18) 

diverges logarithmically in three dimensions as i — > cxd, and a virial expansion 
of transport coefficients does not exist in three dimensions either. Therefore, 
the memory of the initial state is not completely lost. 

The above estimates show that dynamic correlations create long time cor- 
relations so that the initial and kinetic stage are not well separated anymore. 
This observation invalidates the first part of Bogoliubov's assumption about 
the existence of well separated time scales. These long time correlations are 
created by unrestricted free motion, which is possible in isolated groups of 3, 4, 
... particles. However infinitely long free paths between collisions are unphysi- 
cal in a real gas. The free motion of a particle will on average be on the order 
of the mean free path Iq ~ vto ^ l/n. Therefore, the collective effect of the 
other particles, i.e. the surrounding fluid, provides an effective cut off which 
replaces t in (17) and (18) by tg. Consequently, the transport coefficients for 
moderately dense gases exhibit a non- analytic term nlogio ~ M logn| in the 
density, i.e. 

^70 + '71'^ + log ((^ = 3) 



"^-{Vo + V'i^ogn (d = 2) ^'^^ 

where the coefficienls. r][ and 772 have been calculated for hard sphere fluids @ 
and Lorentz|| gases E3. 

■^Thc Lorentz gas is a simple model for diffusion and consists of independent light particles 
moving through a random array of fixed scatterers. 
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The analytical results (19) for transport coefficients can be considered 
as well established for fluids with short range repulsive interactions, and are 
believed to be valid for realistic interparticle potentials asjiiell. The low density 
expansion (19) appears to be consistent with experiment EZI, but so is the virial 
expansion (15). However, up to the present time there are no compelling direct 
experimental observations for the appearance of a logarithmic term log n in 
the density expansion of the transport coefficients. 

6 Necessity of log n singularities? 

The divergence of the ring colhsion integrals, discussed in (17) and (18), holds 
for continuous velocity variables, and so does the logarithmic singularity (19) 
in the low density limit. However in Siistems with discrete vakicity directions, 
such as the Ehrenfest wind-tree modelll3|_ii-dimensional LGAEj, and even one- 
dimensional continuous velocity models E3, the ring collision integrals diverge 
algebraically, instead of logarithmically. This can have very surprising effects. 

Consider, for instance, the self diffusion coefficient. in a one-dimensional 
gas of hard rods, which has been calculated exactlyE3. Calculating D from a 
formal density expansion yields nD — do + di{t)n -I- . . ., where di{t) ~ t for 
large t. By taking the collective effects of the surrounding fluid into account, 
t is effectively replaced by the mean free time to ~ 1/n. The renormalized 
low density behavior of the diffusion coefficient becomes nD = {do + d[) for 
small n, i.e. the ring collisions modify the leading order Boltzmann contribu- 
tion. Exactly the same renormalization of the Boltzmann diffusion coefficient 
through ring cnllisions has been found in stochastic lattice Lorentz gases with 
back scatteringE^. The exact low density diffusion coefficient for these models, 
which is very different from the Boltzmann value, has been calculated by van 
Bcijeren and Ernst EiJ through an exact enumeration method. 

These two examples have been chosen in order to illustrate (i) that the 
behavior of transport coefficients, even at low or moderate densities, is not 
necessarily determined by the dynamics of small isolated groups of 2,3,4,... 
particles, as had been assumed in the Bogoliubov Choh Uhlenbeck theory, 
but that it is strongly influenced by collective effects, and (ii) that the ring 
collisions not necessarily lead to (logarithmic) singularities at small densities. 
It all depends on the structure of the relevant phase space. 

The previous discussion may suggest that in principle the transport prop- 
erties in two- and three-dimensional systems are well understood. We know the 
coefficient of the dominant singularity n''~^logn in d dimensions [d = 2,3). 
We know how to renormalize all divergent terms of the n-expansion into finite 
contributions by resumming the most divergent ring diagrams. It looks as if 
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only technical difficulties are stalling further progress. _ 
However, we do not know the total coefficient of 0{n'^~^). It was PomeauEa, 
who first realized that the total coefficient of the C'(n^)-terni in two dimen- 
sions was divergent as well, and that there was something fundamentally wrong 
with the Navier Stokes equations and transport coefficients in two-dimensional 
fluids, whereas the hydrodynamic equations in three dimensions are still well 
defined to Navier Stokes order. 

To conclude this section, it is worth mentioning that calculation of terms 
beyond the dominant log n-term in three dimensions, for instance the 0{n'^)- 
contribution, is exceedingly difficult, and general results are lacking. In fact, 
no progress along these lines has been made since the early eighties. The sys- 
tematic kinetic theory for fluids has reached an impasse, and the high goals of 
the modern era have unfortunately not been realized. The reader, interested 
in further details, is referred to the review by Dorfman and van Beijerena. Al- 
though twenty years old, this review still gives the best account of the current 
status of the systematic kinetic theory. 



7 Long time tails 

In 1970 Alder and Wainwright^l made their famous discovery of the long time 
tails in the velocity correlation and other current correlation functions in equi- 
librium fluids, although a first indication of the existence of long tirnn^tails in 
three dimensions was already present in a paper of 1966 by GoldmanE^. Com- 
puter simulations on hard sphere systems otj show that the velocity correla- 
tion function (VAF) of a tagged particle, {vx{t)vx(0)), has a positive long-time 
tail, proportional to t^'^/^, which is caused by the following mechanisma. The 
initial momentum of the tagged particle is transferred in part to the surround- 
ing fluid particles. This sets up a flow of vortices around the tagged particle, 
which transfers some of its momentum back to the tagged particle ("kick in 
the back"), thus yielding a positive correlation between v(0) and v(t). This 
correlation extends over hydrodynamic time scales. In LGA similar long time 
tails have been observed with much higher statistical precision and over much 
longer time intervals, using Frenkel's moment propag ation method El. The 
long time tails are a convincing illustration of the important role that com- 
puter simulations have played in unraveling the complex dynamics of classical 
fluids. 

The hydrodynamic picture can be cast into a more quantitative foTHi by 
using the mode coupling theories, developed for critical dynamics E3ij, to 
fluids away from crytical points. As an example consider {vx{t)vx{0)) for 
a tagged particle. To describe the long time behavior on a coarse grained 
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spatial and temporal scale, I decompose the total tagged particle current, 
v(t) ~ J drns{r,t)u{r,t) in terms of the two slow modes: tagged particle 
density ns{r, t) and fluid flow velocity u(r, t). The relevant Fourier modes de- 
cay as ni(k, t) ~ exp{—Dk^t)hs(]s.,0) and u(k, t) ~ exp(— i^fc^i)u(k, 0), where 
D is the coefficient of self diffusion and v — rj/ p the kinematic viscosity. Con- 
sequently, the total current decays as 

vAt) J eM-iD + ^)kH] ^ [A7:{D + v^/^ (20) 

This leads to the quantitative explanation of long time tails, as given by the 
author, Hauge and van Leeuwen E2I. A more fundamental approach through 
resummation of ring collision events was given by Dorfman and Cohen cj. 
The former theories are valid for general densities, the latter are restricted to 
low and moderate densities. I also want to point out that the kinetic theory 
approach of PomeauE3, together with that of Dorfman and Cohen, in which 
the ideas of Pomeau's 1968 articles were incorporated and clarified, might be 
considered as the origin of ring kinetic theory, as it is used at the present times. 

The predicted asymptotic long time tails in the stress correlation function 
in (7) and in the other Green Kubo integrands have never been observed in 
computer simulations in hard core fluids. They are overshadowed by the mo- 
lasses tails, which are present on intermediate time scales, and are related to 
the cage effects ( see section 11). However in LGA with only on-node inter- 
actions the mechanism for "building cages from surrounding hard objects" is 
absent, ancLthe observed long time tails agree reasonably well with theoretical 
predictionsEj. i_. 

In purely diffusive models, such as Lorentz gasesliS, the correlation func- 
tion decays as {vx{t)vAQ)) ^ —t^^^'^^^, i.e. after a long time t the moving 
particle returns to its point of origin on average moving in a direction opposite 
to its initial velocity, instead of parallel as in the fluid case. Hence, the cor- 
relations are negative at large times and weaker than in fluids. In continuous 
Lorentz gases, as well as in LGA the observed long time taijjihfinomena were 
quantitatively explained by means of mode coupling theoryEilcJ as well as by 
ring kinetic theory taO. No new fundamental problems seem to arise here. 

Experimental confirmation of the lopg time tails for air and gaseous argon 
has been obtaincid, by Kim and Matta and for latex spheres in water by 
Ohbayashi et al. Ea using the photon correlation mehtod. Neutron scattering 
experiments on liquid argon and sodium also, offer indirect evidence for the 
existence of long time tails (for references seetl). 
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8 Renormalized fluid dynamic equations 



The most dramatic consequence of the long time tails in fluids is the non- 
existence of the Navier Stokes equations in two dimensions E3. As the trans- 
port coefficients of viscosity, heat conductivity and self-diffusion are expressed 
through Green Kubo formulas as time integral over current correlation func- 
tions, these time integrals are diverging as logi, and transport coefficients do 
not exist in two-dimensional fluids. In purely diffusive systems, like Lorentz 
gases, on the other hand. Pick's law is still valid in two-dimensions, since 
{vx{0)vx{t)) ^ 1/t^ for large and the diffusion coefficient D is finite. 

In three dimensions conventional hydrodynamics involves constitutive re- 
lations where (Fourier transforms of) fluxes are linear in the gradients or in the 
wave number k and the coefficients of proportionality are the Navier Stokes 
transport coefficients. In the same conventional picture one expects that the 
fluxes can be expanded systematically in powers of fc, where the coefficients 
of proportionality in the (!?(fc^)-terms are the Burnett coefficients. This sys- 
tematic expansion in the non-uniformity parameter, fi ~ ^oV, was the basis of 
Bogoliubov's generalization of the Chapman Enskog method in section 3. 

HossseiKr, explicit calculations using mode coupling theory or ring kinetic 
theoryr^rr^ show that the Burnett coefficients in three dimensions are diverging 
logarithmically. This is a consequence of the occurrence of algebraic long 
time tails, which show that the kinetic and hydrodynamic time stages are not 
well separated. This observation invalidates the second part of Bogoliubov's 
assumption about the existence of well separated time scales. 

We observe that there is in fact a close parallel between on the one hand 
the non-existence of the density expansion in (16), which resulted in a non- 
analyticity in the low density behavior of the transport coefficients, and on 
the other hand the non-existence of the /i-expansion, i.e. the non-existence of 
the Navier Stokes equations in two dimensions, and the non-existence of the 
Burnett equations in three dimensions. 

The legitimate question is then: Are there any singularities in the renor- 
malized /i-expansion? If the Navier Stokes equations do not exist in two di- 
mensions, what is the structure of the renormalized hydrodynamic equations 
in two-dimensional fluids? Similarly, what is the structure of the renmrialis 
fluid dynamic equations in three dimensions? Mode coupling theory Ej'Ej'i 
as well as kinetic theory provide partial answers. 

Suppose one calculates the dispersion relation for the relaxation rate z±{k) 
of transverse velocity excitations (shear modes) . Conventional hydrodynamics 
gives z±{k) = —vk^ + ■ ■ where v = rj/p is the kinematic viscosity. The new 
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theory yields: 



( -.P + .[k^f^ + ... id = 3) 

-i^'k^\ogk + ... (d = 2) ^^^^ 

The small-fc singularity in the dispersion relation suggests that the constitutive 
relations between fluxes and-gcadient fields involve transport kernels which are 
nonlocal in space and timeE3cj. There exist also more complex systems with 
intrinsic length scales larger than Iq ( e.g. dense fluids ) where generalized 
hydrodynamics with wavenumber dependent elastic and dissipative reponse 
functions are relevant tl. The same holds for LGAHd. 

The general structure of the renormalized fluid dynamic equations in two 
and three dimensions is still an open problem. A clear resolution of their struc- 
ture is hampered by the apparently complicated non-local structures produced 
by the long range spatial and temporal correlations, as well as by the small- 
ness of the effects that these correlations produce when they are calculated or 
hunted for experimentally. As the effects of long time tails are very small in 
three dimensions, the use of the conventional Navier Stokes equations should 
ofcourse continue in the forseeable future. 



9 Long range correlations in nonequilibrium stationary states 

The first author to analyze already in 1978 long range spatial correlations 
in non-equilibrium fluids in terms of ring collision events was OnukicS. He 
calculated the pair distribution function in a steady laminar flow using the 
Boltzmann-Langevin equation, and showed that the pair correlation function 
has a very long range, exhibiting Coulomb-type behavior, and is proportixjiLal 
to log r for d = 2 and to l/r for d = 3. Kirkpatrick et al. (see review 0c3) 
used in 1982 ring kinetic theory to calculate long rang correlations in fluids 
with an imposed temperature gradient (here the range is even longer than in 
the laminar flow case), and calculated the dynamic strucure factor, which is 
measured in light scattering experiments. 

Such spatial correlations occur in open systems that are kept out of equi- 
librium by reservoirs c3 or by external driving fields IIj. The objects of inter- 
est here are the spatial correlations between fluctuations in the local density 
Sn{r,t) = n{r,t) — (n) or in the local flow velocity du{r,t) = u(r, t) — (u). 
The correlation functions are taken at equal or at unequal times. The spatial 
Fourier transforms of the former are the static structure factors; the spatial 
and temporal Fourier transforms of the latter are the dynamic structure fac- 
tors. The correlations and structure factors related tOj-depsity fluctuations 
are studied in light and neutron scattering experiments c3>El; those related to 
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fluctuations in the flow field are studied in homogeneous turbulence 0. The 
theoretical and experimental results for non-equilibrium stationary states in 
fluids with externalljpmposed gradients in temperature or flow field have been 
recently reviewed inca. Those for driven diffusive systems are discussed in a 
recent monograph by Schmittmann and Ziacj. 

Also here it is instructive to discuss a simple example: I consider the 
density-density correlation function G(r) and corresponding susceptibility or 
structure factor x(k) in the stationary state of a spatially uniform driven sys- 
tem, where the driving field violates the conditions of detailed balance and 
breaks the spatial isotropy, i.e. 

G(r) = - [ dr'{Sn{r' + r)Sn{r')) 



V „ 
1 



dkx(k)exp(ik-r). (22) 



Clearly, for the correlation function G(r) to have an algebraic tail, the sus- 
ceptibility xO^) needs to have a singularity at small wave number k. The 
singularities discussed in section 8 are branch points of the form log k or fc" 
with non-integer a; the same holds in section 7 for the time correlation func- 
tions with r replaced by the time t and k by the frequency lo. However, the 
singularities need not be branch points. 

In the present case of a vector field x(k) the singularity turns out to be a 
discontinuity at k = {k^, ky, ■ ■ ■} — O^-Ihis^is shown by explicit calculations 
based on fluctuating hydrodynamics CjoIIj or ring kinetic theory C3. The 
presence of the driving field makes the susceptibility anisotropic, i.e. 

lim x(k) = Xo(k) (23) 

k — >oc 

depends on the direction k along which the origin is approached. The long 
range part of G{r) is then given by 

G(r) = / dkxoik) expizk • r) = (24) 

where a rescaling of the k-variable shows that the pair correlation function 
decays algebraically, with a coefficient h{r) that only depends onjie direction 
of r. An illustrative landscape plot of G(r) can be found in Refcll. Also, note 
the importance of the anisotropy. If Xo were to be a constant, independent of 
k, then the correlation function above would reduce to G(r) = XoS{r), and is 
strictly short ranged. The algebraic tails ~ l/r'* are intimately related to the 
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long time tails l/f^/"^, discussed above. In fact, they are driven by the same 
slow diffusive modes of the ring kinetic equation, where spatial and temporal 
scales are related by ^ t. 

Also in the context of lattice gas automata the long range spatial corre- 
lations in velocity and density fields have been observed in simulations over 
about 1.5 decades on the spatial scale, and the theoretical calculations ob- 
tained from ring and repeated ring kinetic theory areJji excellent agreement 
with the results from molecular dynamics simulationsllS. 



10 Ring kinetic theory 

The "red thread" connecting all sections in this review are the ring collisions, 
which were born out of the systematic analysis of the triple collision integrals 
of the Choh Uhlenbeck equation. Hence the name "cradle of modern kinetic 
theory" . Ring kinetic theory has been at the basis of most major developments 
in nonequilibrium statistical mechanics over the last 40 years: it explains the 
breakdown of the virial expansion of transpoi± coefficients in continuous flu- 
ids, and their logarithmic density dependenceu; it explains the algebraic long 
time tails of Green Kubo cuxrent correlation functions in continuous fluids □ 
and in lattice gas automata E3; the breakdown of the Navier Stokes equations 
in two dimensions for very long times, as well as the non-analytic dispersion 
relations for sound propagation and for relaxation of hydrodynamic excita- 
tionslj. Moreover, it explains the existence of long range spatial correlations in 
nonequilibrium stationar\!_states, coupled to reservoirs that imiKise shear rates 
or temperature gradientsc3S, or in driven diffusive systems ocJ. 

It seems therefore appropriate to give a siaiple derivation of the ring kinetic 
equation for hard spheres, as given in Refs.L3c3. Under the restriction of low 
densities, the BBGKY hierarchy for hard spheres with diameter a can be 
obtained from (4) and (5) by replacing 9i2 by the binary collision operator: 

Toil2)^S{ri2)a''-^ [ da|vi2 ■ <t|(6^ - 1). (25) 

J vi2 ■a-<0 

The vector cr connects the centers of the two colliding spheres at contact and 
hats denote unit vectors. The integration over the solid angle da is restricted 
to the precoUision hemisphere, vi2 • cr < 0. The operator bcr replaces the 
precoUision velocities {i — 1,2) by postcoUision ones, i.e. 

= &crVi = Vl + (vi2 • a-)a- 

V2 = fo(TV2 = V2 - (Vi2 ■ (26) 
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The Boltzniann equation for a dilute gas of hard spheres follows from the hard 
sphere hierarchy by combining the first equation with the molecular chaos 
assumption (2). As the density increases the ring collisions build up dynamic 
correlations, which invalidate the molecular chaos assumption. 

The method is based on a cluster expansion of the s-particle distribution 
functions, defined recursively as 

/i2 = /1/2 + gi2 

/123 = /1/2/3 + /l 523 + /2513 + /3.gi2 + 3123, (27) 

and so on. Here 1712 accounts for pair correlations, (7123 for triplet correlations, 
etc. The molecular chaos assumption implies gi2 = 0, which is equivalent to 
(2). The basic assumption to obtain the ring kinetic equation is that pair 
correlations are dominant and higher order terms in (27) can be neglected, i.e. 

ffl23 = ffl234 = . . . = 0. 

Setting gi2 = in the first hard sphere hierarchy equation and eliminating 
dfi/dt (z = 1, 2) from the second hierarchy equation with the help of the first 
one, yields the ring kinetic theory of elastic hard spheres: 

{dt + C\l))h^ j dx2To{l2){hh + g^2) 

{9t + /:°(12) - To{l2)-{\ + Vi2) j dx^W?,){\ + Vi^)h}gi2 

= ro(12){/i/2+5i2}. (28) 

Here Vij is a permutation operator that interchanges the particle labels i and 
j. The second equation is the so called repeated ring equation for the pair 
correlation function. If the operator ro(12) on the left hand side of the sec- 
ond equation is deleted, one obtains the simple ring approximation. Formally 
solving this equation for 512 yields an expression in terms of the single-particle 
distribution functions fi (i — 1, 2, 3), and subsequent substitution into the first 
hierarchy equation above yields the generalized Boltzmann equation in ring ap- 
proximation. Essentially all results about singularities and tails, discussed in 
sections 5 to 9, can be obtained from the ring kinetic equation, as given in 
(28). However, the convergent Choh Uhlenbeck correction nrji to the three- 
dimensional viscosity in (15) is not correctly accounted for in this approximate 
equation. 

11 Molasses tails 

At liquid densities the velocity correlation function shows a negative part at 
intermediate times, the so-called cage effected. At the same density and in 
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the same time interval the stress correlation function in (7) appears to have 
in three dimensions a long time tail, that can be fitted to i"'^/^, but with a 
coefficient 400 to 500 times larger then predicted by the mode coupling theory 
of section 7. This so-called molasses tail occurs in dense liquids with short 
range hard core interactions, where the mean free path £o is much smaller 
than the hard core diameter a. 

The conventional long time tail ^ t^'^/^ is controlled by long wavelength 
modes, and the intermediate time molasses tail by slow modes of intermediate 
wavelength around ka ~ 27r, where the structure factor has its maximum. At 
these wavevectors the dynamics of fluid excitations in fact controlled by a 
single soft (density) mode, also called cage diffusionu. At the high densities 
relevant for liquids, the Fourier modes of momentum and energy density with 
wave vectors near 27r/tT relax very rapidly through the mechanism of instanta- 
neous collisional transfer over distances of order cr, whereas kinetic transport, 
which is operating on scales of £o, is much slower. Therefore the density, which 
can only relax through kinetic transport, is the only slow mode near ka ~ 27r. 
It requires rearrangments of particle configurations over distances of order a 
(see reviewlj), i_e. structural relaxation. The detailed explanation was given 
by KirkpatrickEd and van BeijerenEJ, who extended the mode coupling theory 
of section 7 by including wavenumbers upto k ~ 27r/(T, i.e. by including the 
slow cage diffusion mode. This important fundamental result marks perhaps 
the end of the modern era. 



12 Post-modern kinetic theory 

By 1985 the results of the modern era: the log n-dependence of transport 
coefficients, the long spatial and temporal tails in correlation functions, and the 
fact that conventional hydrodynamics breaks down, at least in principle, were 
generally accepted. However, the experimental consequences of the exciting 
developments were minute and difficult to observe experimentally. 

Where to go from here? The tendency was to combine intuitive classical 
concepts with the modern concepts of dynamical correlations, originating from 
mode coupling and ring kinetic theory, and to apply these to, more complex 
problems, such_as nonlinear hydrodynamics and porous mediaEI, dense fluidsB, 
granular flows e3, mathematical modeling of the Boltzmann equation and ring 
kinetic equations Ej, and even to the calculation of Lyapunov exponents in 
non-equilibrium fluids O. Simultaneously one sees the development of new 
methods for numerical|-simulation of non-equilibrium fluids, such as Lattice 
Gas Automata (LGA) eL the Lattice Boltzmann Equation (LBE) method E3, 
and Bird- type methods E3, which are stochastic simulation methods for solving 
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nonlinear kinetic equations. 

Many of these current developments, such as the LBE method, the Bird- 
type methods, mathematical modeling of kinetic equations, and Lyapunov ex- 
ponents, have no direct links with the modern concepts of dynamical correla- 
tions and with the Choh Uhlenbeck equation, and will not be discussed any 
further. 

An outburst of renewed interest in kinetic theory occurred in 1986, when 
Frisch, Hasslacher and PomeauEl proposed Lattice Gas Automata (LGA), as 
models for non-equilibrium fluids. This event marks in my view the beginning 
of the Post-modern Era. The LGA, already defined in section 2, are fully 
discretized statistical mechanical models of A'^-particle systems, which do not 
involve any new fundamental concepts. The LGA were specially designed for 
large scale and long time computer simulations of complex hydrodynamic flow 
problems. As the dynamic description is based on Boolean variables, there are 
no round off errors. The kinetic theory for these models closely parallels that 
for continuous fluids, as I have been indicating in all previous sections. 

The most important fundamental result from LGA is the excellent con- 
firmation of the theoretical predictions for the long time tails in the velocity 
correlation function, calculated from mode coupling and ring kinetic theory. 
Here the tails have been meiasured with very high statistical accuracy over two 
decades of mean free times cj. 

Moreover, the study of LGA has provided many interesting new results on 
the dynamics of phase separation and pattern formation, as well as on fiows 
through porous media (see review by Rothman and ZaleskiEl). Unfortunately, 
LGA turned out to be not very suitable to study high Reynolds number flows, 
because the viscosity of LGA is too large. 

A phenomenological equation linked to the modern era, and widely applied 
in the post-modern era, is the nonlinear revised Enskog equation for hard 
spheres, which is a MaiJiovian approximation, constructed in 1972 by van 
Beijeren and the author EH as an exact result for short limes. The linearized 
version of this theory has been derived by many authorsl£3. These theories are 
collectively referred to as the Revised Enskog Theory (RET) with a collision 
operator taking non-local short range static correlations into account. 

The theoretical and experimental research by de Schepper and coUabora- 
torslj on dynamical structure functions of dense liquids with short range hard 
core interactions has shown that the RET gives a reasonably accurate descrip- 
tion of the elementary excitations with wavelengths in the range £o < A < ct. 
A Fourier mode analysis Ej of the revised Enskog equation shows that this 
equation has only a single slow mode at wavenumbers near fccr ~ 27r, the cage 
diffusion mode. Although the RET involves only binary collision dynamics, it 
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mimics the effects of cage diffusion, described in section 11, reasonably well. 
The classical Enskog equation in (3) is missing this effect. 

Further interesting extensions of the RET to non-equilibrium solids have 
been given by Kirkpatrick et alL3, who derived the macroscopic equations of 
motion for the elastic solid from this equation. RET yields in principle a mi- 
croscopic theory for elastic constants and dissipative coefficients, as well as 
a dynamic theory for freezing. In the fluid phase there are five slow hydro- 
dynamic modes; in the solid phase there are eight soft modes with vanishing 
dispersion relations as k 0, i.e. five conserved densities and three displace- 
ments fields (order parameter fields). There are six propagating modes, a ther- 
mal diffusion mode and a vacancy diffusion mode. Vacancy diffusion, which 
involves structural relaxation, occurs on time scales much greater than those 
occuring in sound phenomena and thermal diffusion, which are controlled by 
collisional transfer of momentum and energy. This suggests that cage diffusion 
in the liquid phase and vacancy diffusion in the solid phase may be related to 
the same slow microscopic excitation. Moreovp:. the RET has also been ap- 
plied by KirkpatrickEil, Kawasaki e3 and DuftyEa to study the slow structural 
relaxation in glasses. 

The ring equation of section 10 has still a very complex structure. Dufty, 
Lee and BreyE3 have proposed a mathematical model with a simpler structure 
than (28). It would be very interesting to obtain from this simpler model more 
explicit analytic results for the pair correlation function injion-equilibrium 
fluids. Finally, a very recent publication by van Noije et al.E3 on long range 
correlations in rapid granular flows shows that ring kinetic theory or fluctuating 
hydrodynamics has become the standard tool to analyze long range correlations 
in complex and simple fluids out of equilibrium. 

The central theme of this review has been the importance in non-equilib- 
rium statistical mechanics of the ring collisions, which were born out of a 
detailed analysis of the Choh Uhlenbeck equation. 
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